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Abstract 

Decay rates of the B-meson are studied through charmless and strang- 
less transitions into tt, p , u and 7 systems. The important features of 
these modes are their clean signatures. The CLEO II collaboration has 
recently reported the value BR(B° -> vr+vr^) = ± 0.2 x 10~ 5 . 

This value test our approach to nonleptonic B-decays within the stan- 
dard model (SM). Since B — > 77 is far beyond our experimental reach, 
we believe that the correct determination of the order of magnitude 
~ 10 for BR(B — > 77) provides the most reliable value needed at 
this moment. The most recent experimental report by the CLEO 
collaboration on BR(B -> K*i) = (4.5 ± 1.5 ± 0.9) x 10~ 5 rep- 
resents a confirmation of our SM prediction for the B — > K*j de- 
cay. This experimental result, despite of its great importance for the 
SM and the physics beyond the SM, enables us to predict the value 
BR {B -> py) = (3.5 ± 3.3) x lO" 6 . 



*e-mail: josipt@rex.irb.hr 



1 Introduction 



Hadronic rare B decays have recently been the subject of both theoretical 
and experimental interest. The estimates for these processes involve the 
matrix elements of four quark operators, and these did not receive much 
attention in the past. In this paper we propose to study a number of B 
decays that follow from b — > u,d and have clean signatures. The calculation 
proceedes in two steps. The effective short- distance interaction consists of 
two contributions. One contribution comes from the QCD-corrected tree- 
level w-exchange diagram [1], with the b — > u vertex. The other contribution 
comes from the gluon exchanged with b — > d transitions [2]. The second 
step is to use the factorization approximation to derive the hadronic matrix 
elements by saturating H%ff with the vacuum state in all possible ways. 
The resulting matrix elements involve quark bilinears between one meson 
and the vacuum and between two meson states. These matrix elements are 
estimated using relativistic quark model wave functions. Such a technique 
has been extensively used for B and D nonleptonic decays by Bauer et al [3] 
and the results are consistent with experiment. These methods for b — > s 
nonleptonic modes were employed [4] in our previous study. 

The short- and long-distance contribution to B — > p7 and the determi- 
nation of the order of magnitude for B — > 77 are described in the last two 
sections of this paper. 

2 Effective Hamiltonian 

The effective weak interaction Hamiltonian (local four-quark operator) de- 
scribing the b — > u transitions for charmless and strangeless B-meson decays 
at the tree level is 

H% f = V2G F (c-0- + c + 0+) + h.c, (1) 

Ot = ^A^Xi^A^X- (2) 

Here c± are the QCD coefficients and (i, j) are the color indices. The param- 
eters used to evaluate the c± coefficients are 

m b = A.9GeV, A QCD = 200MeV, a s (m b ) = 0.256, (3) 
c_ = 1.34, c + = 0.86. (4) 
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Charmless and strangless B decays can also arise through a one-loop process 
involving gluon exchange. The relevant Hamiltonian is [2] 

+ tiA^L-^K^ti), (5) 



o> s GfGi g 1 

Here q runs over all quark species, although only u and d are relevant to our 
discussion. We have used a s (mb) = 0.256 as in (J3J), and G\ = —6.52Vb c V* c 
corresponding to m t = 150GeV. 



3 Factorization approximation 

From experience we know that nonleptonic decays are extremely difficult to 
handle. For example, the AI = 1/2 rule in K — > ttti decays has not yet 
been understood in a satisfactory way. Enormous theoretical machinery has 
been applied to K — > nn decays producing only up to 50% agreement with 
experiment. For energetic decays of heavy mesons (D, B), the situation is 
somewhat simpler. For these decays, the direct generation of a final meson 
by a quark current is (probably) a good approximation. 

According to the current-field identities, the currents are proportional to 
interpolating stable or quasi-stable hadron fields. The approximation now 
consists only in taking the asymptotic part of the full hadron field, i.e. its 
"in" or "out" field. Then the weak amplitude factorizes and is fully deter- 
mined by the matrix elements of another current between the two remaining 
hadron states. For that reason, we call this approximation the factorization 
approximation. Note that in replacing the interacting fields by the asymp- 
totic fields, we have neglected any initial or final-state interaction of the 
corresponding particles. For B decays, this can be justified by the simple 
energy argument that one very heavy object decays into two light but very 
energetic objects whose interactions might be safely neglected. Also, dia- 
grams in which a quark pair is created from the vacuum will have small 
amplitudes because these quarks have to combine with fast quarks to form 
the final meson. Note also that the 1/N C -expansion argument provides a 
theoretical justification for the factorization approximation, since it follows 
the leading order in the 1/N C expansion. Here N c is the number of colors. 
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Each of the B-decay mode might receive three different contributions. As 
an example we give one amplitude obtained from H™ff '■ 

A(B + (p)^7r + (k)n\q)) = L(n°)(n + |%(l - l5 )d\ j B + )(tt° \u^ l5 u\ 0) 

+ L(7r+)(7r+ \u^j 5 d\ 0)(tt° |& T "(1 - j 5 )u\ B + ) 
+ L(B + )(n + n° \u lp (l - l5 )d\ 0)(0 \h^\ B + ). 

(7) 

The coefficients L(tt°), L(tt + ) and L(B + ) contain the coupling constants, 
color factors, flavor symmetry factors, i.e. flavor counting factors and fac- 
tors resulting from the Fierz transformation of the operators in Eqs.Q and 
(JH). The coefficients L(tt°) and L(n + ) correspond to the quark decay dia- 
gram, whereas the L(B + ) corresponds to the so-called annihilation diagrams. 
These factors are different for each decay mode, as indicated by the depen- 
dence on the final-state meson. To obtain the amplitudes for other decay 
modes, one has to replace the final-state particles with the particles relevant 
to that particular mode. 

The QCD coefficients appear in two different combinations in the ampli- 
tudes of various decay modes [5]: 




The factor 1/N C arises from the color mismatch in forming color singlets after 
Fierz transformation. 

We proceed with the definitions of the coupling constants and the Lorentz 
decomposition of the typical hadronic matrix elements: 

(0 \h^ 5 u\ B+(p)) = ip„f B , f B = 1.5f w , U = 130MeV, (10) 
(p°(k) |m7 m m| 0) = ig p oe p {k), (u(k) |«7 M w| 0) = ig^e^k), g p o = (11) 

(n + (k) \h,d\ B + {p)) = (p + fcW (+ V) + <lJ H (q 2 ), q = P~k (12) 

(p + (k) \h,(l - 7sH B + (p)) = is, uX y{k){p + k)\p - kyv{ q 2 ) 
+ s,(k)[(m 2 B -m 2 p )A 1 (q 2 )-(E-q)(p + k),A 2 (q 2 )} (13) 

+ q^(^-)(m B + m p ) [2m p A (q 2 ) - (m B - m p ) (A^q 2 ) - A 2 (q 2 ))] , 
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(n+(k)p (q) \u^d\ 0) = (p°(q) \u^d\ 7r~(— &)}. (14) 

Any hadronic matrix element needed to evaluate the branching ratios of 
the decays can easily be obtained from the above definitions. We assume 
that the momentum dependence of the form-factors f + {q 2 ), V(q 2 ), etc. from 
Eqs. (jllHr3j) . is well described by single poles with masses of excited b-quark 
meson states (1~,0~, ...) close to m B , i.e. we can use only the nearest pole 
dominance. The form factors at zero momentum transfer were calculated us- 
ing the relativistic oscillator wave functions [3]. Now we give a few examples 
of the structure of the form-factors: 



f(+) (rr,2\ - fe^j 



(15) 



VK) = V f } zrr , (16) 



(m B + nip) 



< B H) = : ^ =rr . ("> 

(m B + m p ) (1 



< B K) 



-A? (0) 



(m B - m p ) (l - £5 



i+ 



< B «) - . / f l°' -'V (19) 

{m B + m p ) \}-^) 

Here we have used m -(b) ~ m +(b) ~ m>i-(b) ~ m!+(6) »s mg. The 
couplings g p , g u are determined from the p, u — > e + e~ experimental rates: 

* = w r(p ~" e+e_) = °- 0141Gw4 ' ( 2 °) 



27?77 3 

^ = j^^^-e + e-) = 0M28GeV\ (21) 
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For the ratio of the Kobayashi-Maskawa (K-M) matrix elements we use 
the central value |Kb/Hc| = 0.1C, where < ( < 1. Some examples of the 
branching ratios that arise from the tree-level Hamiltonian are 



BR(B° -> 7T+7T-) = 2 

BR (B° - p°cu) 
BR (B^ 



m B 

2 



V, 



ab 



v bc 



\Ci\ 



m B 



V u 



P + ") = ( 



m B 



2 



|C|»(l + ^)V|Ajf. 



where 



Aab — 



m. 



u6 



6c 



d - — c 2 

9ffi 



|F| A 3/2 



2 \ 2 

mi 



\m\m\ 
mi 



(22) 

(23) 
(24) 

(25) 



|F| = 2V 2 (m 2 p ) 

The results for all the modes using the tree-level Hamiltonian are given in 
Table 1. The values of the branching ratios presented in Table 1 have been 
obtained by applying the N c — > oo limit, i.e. for C\ = 1.1 and C 2 = —0.24. 
We find that the penguin contributions to these decays are generally small. 
We therefore estimate the contribution to a few of the more interesting modes 
given in Table 1. Finally note that the branching ratios BR are evaluated 
by normalizing the partial widths V to the B-meson lifetime, i.e. to the total 
B-meson decay width (for details, see Refs. [8,9]): 

„ / „\ 1 G%m% \V r h\ , _x 

r-W = - * ^ V 2 | 3 • < 27 > 
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TABLE 1 

The branching ratios of various charmless and strangless nonleptonic IB- 
decay modes. 



Branching 




QCD coeff. 


Contribution 


Contribution 


ARGUS/CLEO 


Ratio 


Ref.[3] 

[io- 5 ]C 2 


[for the tree 


from H^j: 


from H^ff 


Exp. limit [6] at 
90 % CL [IO" 5 ] 


Mode 


level graphs] 


[io- 5 ]C 2 


[io- 5 ] 


B+ -> 7T+7T 


0.6 


C\ — C 2 


0.49 


0.00 


50/260 


7T+p° 


0.6 


c 2 


0.18 




19/17 


7Y + U 




c 2 


0.18 






n°p+ 


0.6 


Ci 


1.74 






P + P° 


1.4 


Cl — C 2 


1.15 






p + u) 




Eq.(24) 


1.19 






B° -> 7T + 7T- 


2.1 




1.56 


3 x 10~ 4 


1.3^g ±0.2 [7] 


7T 7T 




c 2 


0.03 


2.5 x 10~ 4 




7T + p~ 


5.6 




4.35 










Ci 


1.07 






7T p° 


0.1 


c 2 


0.10 


1.8 x 10~ 3 


43 






c 2 


0.09 






p + p~ 


4.5 


Ci 


3.60 


4 x 10~ 3 


420 


p°p° 


0.1 


c 2 


0.07 






p u 




c 2 


0.14 






Ijjljj 




c 2 


0.06 







4 Short and long distance contributions to 

B -> pj 

In April this year the CLEO collaboration announced the experimental dis- 
covery of electromagnetic-penguin B decays B — > K*^. The information 
from Wilson Lab on 13-APR-1993 was the following [10] : 

Using a data sample of 2.8 x 10 6 B-meson decays collected by the CLEO II 
detector operating at the Cornell Electron Storage Ring (CESR) it was possi- 
ble to observe the rare b-quark decay b — > s'j in the modes -B — > K*° (892)7 
and B~ — > 7^*^(892)7 (the charge conjugation was implied throughout). 
Details are given below in Table 2. This is the first unambiguous evidence 
for penguin-type diagrams in weak decays. 
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TABLE 2 



The CLEO collaboration experimental discovery of electromagnetic- 
penguin B-decays [10]: B — > K*j. 



Mode 


K* u 7 


K* 


7 


K* decay mode 


K 7T+ 






Signal events 


8 


2 


3 


Sideband events 


41 


2 


10 


Sideband scale factor**) 


1.0/37.6 


1.0/40. 


1.0/12. 


Extrapolated background 


1.1 ±0.2 


0.050 ± 0.035 


0.83 ±0.26 


Binomial probability of 
background fluctuation 


3.7 x 10~ 5 


0.35% 


7.30% 


Efficiency 


(11.9 ± 1.8) % 


2.0% 


3.1% 


Additional BB backgnd ( ** } 


0.30±0.15 


0.01 


0.10 


Branching ratio [10 -5 ] 


4.0 ±1.7 ± 0.8 


5.7 ± 3.1 ± 1.1 


Combined BR [B -> (K M + K*~)-f] 


(4.5 ± 1.5 ±0.9) x 10~ 5 



(*)This is the number of background events in the signal region 
divided by the number of background events in the sideband 
region from off BB resonance data. 

(**) This is the number of background events expected from BB 
events from the Monte Carlo simulation. Note that the background 
is dominated by non-BB processes. 

The branching fractions were calculated by summing signal events, sub- 
tracting summed background events and dividing by summed efficiencies. 
In particular, the combined K*° and K*~ number was NOT obtained by 
averaging the individual results, weighted by the statistical (or any other) 
error. The combined number was obtained by summing all signal events, 
subtracting all backgrounds and dividing by the sum of all efficiencies. 

The short distance (SD) contribution to B — > K* r y was proposed and 
obtained in several papers [11,12]. The branching ratio obtained is of the 
form 



BR SD (B -> iT 7 ) = 



a / rrib 
2tt \m B 



v t * s v t 



ib 



V, 



be 



jf K \0)\ 2 F?(m t ) . (28) 



Here F$ s (mt 



is the one-loop flavor- changing b — > 57 vertex with QCD cor- 
rections, and jf K (0) is the form-factor of the operator sa^q^bu between 
the spin-one kaon resonance and the B meson, responsible for the b — > S7 
transition at q 2 = o, i.e. for the real photon. Details of the calculations are 
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given in Ref. [8]. The BR SD is obtained by normalizing the partial width 
T SD (B —>■ K*j) to the B- meson lifetime (see Eq. (I28J1 ). For the top-quark 
mass m t = 150 GeV, we have found [8]: 



BR SD (B -> K*i) = 2.6 x 10" 



(29) 



which is in excellent agreement with the recent theoretical estimate [13] and 
is within the error of the measured branching ratio published recently by the 
CLEO collaboration [10]. 

The short distance contribution to B — > /ry is dominated by the flavor- 
changing vertex b — > d'-f , which proceeds in one loop through the exchange 
of u,c and t quarks and the W-boson. For the emission of the real photon 
the only contributing operator da^ u q v b R gives the following branching ratio: 



BR SD (B -> pj) 



at_ I m b y 
2tt \m B J 



VuV tb 



Vi 



be 



fi"(0) 



{m t ) 



(30) 



Here F 2 M (mt 



is the one-loop flavor-changing b — > d'j vertex with QCD cor- 
rections, and the fi P (0) is the form-factor of the operator d<T ia ,q u bn between 
the p-meson and the B-meson. The BR SD is obtained by normalizing the 
partial width T SD (B -> fry) to the B -meson lifetime. 

Careful analysis of the procedure for the one-loop QCD calculations of 
the F2 S (m t ) coefficient [11,12] shows that up to the mass symmetry-breaking 
effects, one can write 

F b 2 s {m t ) = F h 2 d (m t ). (31) 

The calculations of the matrix element (K* + (k) \~s<y fJ-u q u bR\ B + (p)), in the 
Relativistic Constituent Quark Model (RCQM) [8] gives approximately the 
same value as the calculation of this matrix element [13] in an effective chiral 
theory for mesons using flavor and spin symmetries of the Heavy Quark 
Effective Theory [HQET]. Detailed analysis of both calculations shows that 
up to the quark mass difference m s — m^, (which is negligible compared to 
the b-quark mass), we have 



(K* + (k) \sa^q% R \ B + {p)) = (p + {k) \da, v q v b R \ B + {p)). 
Equation (|3*2"|) in the RCQM then gives 

/f**(0) = /?'(()), /^(0) = /?(0), 



(32) 



From Eqs. 1)311 )1 . ipTjl and it is easy to find that 



BR SD {B — > fry) 



td 



ts 



BR SD (B -> K*i). 



(33) 
(34) 
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The ratio of the K-M matrix elements [9] [V^/V^l = 0.10 to 0.33 and the 
last line in Table 2 gives the following range of values for BR SD (B — > p 7 ): 

BR SD (B^p-i) = (4.5 to 50) x 10" 7 . (35) 

It is appropriate here to comment on the long-distance (LD) contribution 
to B — > fry. The long-distance contribution to B — > K*^ is discussed in detail 
in [14]. Despite of the conclusion in Ref. [14] that the rate for B — > K*ip 
cannot be used to give a unique value for B — > K*^y, it is necessary to evaluate 
BR LD (B — > p 7 ) by applying exactly the same procedure, as in Ref. [14], as 
a check of consistency. We use the vector-meson dominance which leads to 
the vector-meson - photon conversion mechanism. This mechanism is known 
as a long-distance effect. 

Owing to the gauge-condition and current-conservation requirements, for 
the real outgoing photon in the B — > p 7 mode we have that [8] A = and 
A2 = A\ = A. This gives 



BR ld f B o ^ p o 7 ) 




BR LD (B + -> p+ 7 ) 



™9 P + 
m 2 p m 2 B 



c x 1 



ub 



V bc 
2 



i-4)V(o)+^(o)) 



9 p o 



c 2 1 1 + aW > 



(37) 



In the N c — > 00 limit, i.e. for C\ = 1.1 and C2 = —0.24 and for | V u b/Vb c \ 
0.1, we have 



BR LD (B° -> p° 7 ) = BR LD (B° -> cj 7 ) = 0.02 x 10" 



(38) 



BR LD (B + -> p+ 7 ) = 1.4 x 10~ 8 . (39) 

Assuming that the recently quoted result [15] \V u b/Vb c \ — 0.08 is of reasonable 
accuracy, the above BR LD {B — > p 7 ) should become smaller by at least a 
factor of two. This smaller branching ratio, together with Eq. (j35|) . gives 

BR SD (B -> pj) ~ 60 x BR LD (B -> p 7 ). (40) 
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This is similar to our previous result [14] obtained with the top quark mass 
m t = 150GeV for B -> K* T - 

BR SD (B -> K*i) ~ 50 x BR LD (B -> fr 7 ), (41) 

In the best case, the LD contributions to the total branching ratios represent 
~ 2 % corrections for both B — ► fT*7 and 5 — > P7. 

The similarity of the results in Eqs. (jjOj) and (j4*Tj) shows that the pro- 
cedure used in this paper to obtain BR LD (B —>■ fry) is consistent with the 
procedure in Ref. [14]. This means that all conclusions drawn in [14] are di- 
rectly applicable to this work and that the rates for B — > pp,pu,u>ui cannot 
be used to give a unique value for B — > fry, i.e. Eq.fSHJ) fully represents our 
prediction for the decay rate: 

2 x 10" 7 < BR(B -> fry) < 7 x 10~ 6 . (42) 

Here we have also take into account the errors in the measurements [10] of 
B -> K*-f from Table 2. 



5 B decay into the two-photon final state 

In this case, we should only estimate the order of magnitude for BR(B — > 
77). 

First, applying the vector- meson - photon conversion mechanism to the 
neutral B decay modes (p°p°, p u, uu), we may evaluate BR LD (B — > 77) as 



B R LD (B° - 77) 



m 2 B mA 



V, 



be 



(43) 

For C2 = —0.24 (in the iV c — > 00 limit) and \ V u b/Vb c \ = 0.1, we have 

(5° -> 77) = 10" 12 . (44) 

Next, comparing the experimentally measured ratio of widths for if s de- 
cays [9] with the ratio of our theoretical estimates from Eqs. (J38l44|) and Table 
1, we find that 



[T(K S -> 77)/r(K s -> 7r+7r-)] exp = 3.5 x 10- 6 , 



(45) 



[r(B° -> 77)/r(B° -> 7r+7r-)] ift 6.4 x 10- 6 . (46) 
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Since the two rates are in reasonable agreement, it is justified to use the 
conversion mechanism to estimate the order of magnitude for the B° — » 77 
rate to be 

BR (B° -> 77) ~ 1(T 10 . (47) 



6 Discussion and conclusion 

One can see from Table 1 that our values are considerably lower than the 
present limits obtained at 90 % CL. On the other hand, our results are 
in reasonable agreement with the predictions of Ref. 3 summarized in the 
second column in Table 1. Note also that the rates involving the (helicity- 
suppressed) combination (c_ — 2c + ) 2 (or (c_ — c + ) 2 in the N c — > 00 limit) are 
sensitive to the values of the QCD coefficients. Some experimental limits [6,7] 
are also presented in Table 1. 

The CLEO II (Ref. [7]) collaboration has recently reported the BR(B° —* 
7T + 7r~) < 3.0 x 10~ 5 at 90 % CL. This is a significant improvement upon the 
old results (CLEO 1.5) 

BR(B° — > 7r + 7r~) < 9.0 x 10 -5 and is welcome because it is a test of our 
approach within the standard model. 

Before discussing B — > p'-f, let as briefly comment on B — > 77. Since 
the B — > 77 decay is far beyond our experimental reach, we believe that the 
correct determination of the order of magnitude ~ 10~ 10 for BR(B — > 77) 
provides the most suitable value needed at this moment. 

The most recent CLEO report [10] on BR(B -> K*^) = (4.5 ± 1.5 ± 
0.9) x 10~ 5 represents a confirmation of our SM prediction for the B — > K*y 
decay given six years ago [8]. This experimental result, despite of its great 
importance for the SM and the physics beyond the SM, enables us to predict 
the following value for the B — > pry decay within the SM: 

BR (B -> pry) = (3.5 ± 3.3) x 10" 6 . (48) 

The range of the order of magnitude in Eqs. (j42j) and/or ()48|) is due to the 
range of values in the |Vtd/Vt s | ratio [9]. 

Assuming that some other experiment might determine the ratio | V t d/V ts \ '- 
0.1, then at least 10 8 B°B° pairs have to be produced at CESR and collected 
by the CLEO detector to discover the electromagnetic penguin B-decays 
B — > p7. It is clear that 10 8 B-meson decays cannot be collected in the 
near future. However, it is encouraging that Cornell has an asymmetric B 
Factory [16] proposed further CESR upgrade. 

On the other hand, if \V td /V ts \ turns out to be = 0.3, then only 10 7 B°B° 
pairs have to be collected by the CLEO detector to discover the B — > pry 



11 



decay. This scenario is much more promising because Cornell is currently 
upgrading the luminosity of CESR by at least a factor of five [17]. 

In the case of known top-quark mass, the measurement of BR(B — > pj) is 
an excellent tool for fixing the weak flavor-mixing parameters more precisely, 
such as the (V^/V^l ratio. 

As is well known, radiative B-meson decays are of importance as a test of 
the SM electroweak theory in one loop, because they verify the gauge struc- 
ture of the theory. In addition, these decays might contribute to the discovery 
of new physics beyond the SM. Finally, note that the relative smallness of the 
branching ratio for B — > ^7 makes it more suitable for possible discovery of 
new physics than might be expected from the recently measured branching 
ratio for B — > K*j. 
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